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ABSTRACT 

The interaction between an accreting satellite and the Andromeda galaxy (M31) has been studied 
analytically and numerically, using a high-resolution TV-body simulation with 4 x 10^ particles. For 
the first time, we show the self-gravitating response of the disk, the bulge, and the dark matter halo 
of M31 to an accreting satellite. We reproduce the Stream and the shells at the East and West side 
of M31, by following the evolution of the collision 4 Gyr into the future, and show that recently 
discovered diffuse arcs on the south side of the minor axis are the remnants of a similar collision that 
occurred 3-4 Gyr earlier than the Stream event. The present day integrity of the M31 disk constrains 
the satellite mass to be M < 5 x IO^Mq. The stars that were originally in the center of the satellite 
are now in the east shell. Therefore, observations in this region might reveal additional clues about 
the nature of satellite, such as the central core and any metallicity gradient. 

Subject headings: galaxies: interactions — galaxies: kinematics and dynamics — galaxies: individual 
(M31) 



1. INTRODUCTION 

Large spiral galaxies like the Andromeda galaxy (M31) 
are believed to have formed in part from the merger of 
many less massive galaxies, and their halos are modeled 
as the debris of accreted dwarfs (e.g. Bullock & John- 
ston 2005). The giant stream of Ibata et al. (2001) and 
the starcount maps of Ferguson et al. (2002) establish 
the complex merger origin of the M31 halo population; 
evidence for coherent structures now extends to « 100 
kpc (Ibata et al. 2007). Further support comes from the 
spectroscopic side e.g. the disk population of Ibata et al. 
(2005). Deep HST/ACS photometry shows that three 
deep fields in the M31 halo exhibit unequivocal evidence 
for intermediate age populations (Brown et al. 2006). 
Starcount maps and radial velocities of red giant stars 
near M31 exhibit a giant stellar stream to the south of 
this galaxy, as well as giant stellar shells to the east and 
the west of M31's center (Ibata et al. 2001; Ferguson et 
al. 2002; McConnachie et al. 2003; Ibata et al. 2004, 
2005; Guhathakurta et al. 2006). Recently Ibata et al. 
(2007) have carried out a deep, wide field photometric 
survey of the minor axis toward M33. They find a series 
of arc-like structures parallel to the major axis of M31 
toward M33, and spaced by k, 10 kpc. 

A''-body simulations of the interaction between the pro- 
genitor of the giant stream and M31 suggest that the 
stream and shells are the tidal debris formed in the last 
pericentric passage of a satellite on a radial orbit (Ibata 
et al. 2004; Font et al. 2006; Geehan et al. 2006; Fardal 
et al. 2006, 2007). However, these models have always 
assumed a fixed potential to represent the influence of 
M31, and no models consider the possible effect of a live 
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(A^-body) disk. The discovery of the widespread inter- 
mediate age halo population has raised the question of 
whether these stars belong to the infalling satellite or are 
instead disk stars that were ejected in the stream colli- 
sion or an earlier interaction. Present constraints on the 
progenitor mass are weak; consequently we look to the 
present day thickness of the disk to limit the progenitor 
mass. In addition, the newly discovered large scale arcs 
(Ibata et al. 2007) have not yet been modeled. These 
issues motivate us to explore the effect of the accreting 
satellite using the first self-consistent, Af-body model of 
M31 that has a disk, bulge, and dark matter halo. 

2. DISK THICKNESS CONSTRAINS THE SATELLITE MASS 

It has long been known that infalling satellites are 
likely the major source of the observed thickness of disks, 
and that the integrity of disks constrain the masses of 
such satellites (Toth & Ostriker 1992; Quinn, Hern- 
quist & FuUagar 1993; Velazquez & White 1999; Font 
et al. 2001; Hayashi & Chiba 2006; Gauthier, Dubinski 
& Widrow 2006). 

We estimate the disk heating by dynamical friction, 
which is caused by the scattering of the disk stars into 
an overdense wake that trails the orbiting body, exerting 
a force opposite to the orbital motion. The energy in- 
put into disk stars through the interaction between disk 
stars and a satellite is equal to the orbital energy loss 
/S.E of the satellite. If we assume that the vertical veloc- 
ity dispersion of the disk stars CTz is constant, the change 
AcTj. is simply denoted by Acr^ = 2A£'/Ms, where Ms is 
the mass of the satellite, and AE is given by the inte- 
gral of the frictional force over an orbit of the satellite 
within the disk. Using the (Chandrasekhar 1943) ex- 
pression for dynamical friction, for a thin disk, l^E is 
l\E ~ 47r G^Afs^E/Ug , where G is the gravitational con- 
stant, £ is the surface mass density of a disk and Vg is 
the satellite velocity. 
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For an axisymmctric thin disk, the equation of motion 
in the vertical direction and the Poisson equation are 
simply expressed as 9z(pdO'z) + Pddz'^ ~ and d'^'^ — 
AirCpd, where pd is the mass density of the disk and $ 
is the gravitational potential, respectively. Using these 
equations and assuming the pd oc exp(— z/zd), where Zd 
is the scale height of the disk, we obtain the relation 
between a^ and Zd'- <t^ = 27rG'Szd- Inserting the equation 
of AcTg into the change of the vertical velocity dispersion 
of the disk as Act^ = 27rGEAzd, we get the relation as 
Ms = [MdVsAzd/(4G)]^/^. Using a set of parameters 
given by Model A of Widrow, Perrett & Suyu (2003), 
finally, we obtain the critical satellite mass. 



Mc = 5.2 X lO^Mo 
Md 



7 X 1010 Mp 
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150 km s 
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0.3 kpc 
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Consequently, the dynamical mass of the satellite 
should be smaller than this critical mass, since the disk 
thickness must agree with the observed thickness of M31 
after the interaction of the satellite. This yields an up- 
per limit for the total progenitor mass. Adopting for 
the stream [Fe/H]> -1 (Koch et al. 2007) and us- 
ing the mass-metallicity relation of Dekel & Woo (2003) 
gives a lower mass limit for a progenitor stellar mass 
« 5 X 10^ Mq, which agrees with the predictions of Font 
et al. (2006, 2007). We offer only order of magnitude 
constraints on these mass bounds, which must be refined 
further with better observations and fully nonlinear nu- 
merical simulations. 

3. NUMERICAL MODELING OF THE INTERACTION 
BETWEEN THE SATELLITE AND M31 

In this section, we demonstrate an iV-body simula- 
tion of the interaction between an accreting satellite 
and M31. We assume that the total mass of the disk 
is 7.0 X 10i°Mo. The density of the disk falls off ap- 
proximately as the exponential with the scale length of 

5.4 kpc in the radial direction, and as sech^ with the 
scale height of 0.3 kpc in the vertical direction. The 
bulge is a King (1966) model with the total mass of 

2.5 X IO^^Mq. The dark matter halo is taken to be a 
lowered Evans model (Kuijken & Dubinski 1995) with 
the total mass of 3.2 x IO^^Mq and the tidal radius of 
80 kpc. This set of parameters corresponds to Model 
A of Widrow, Perrett & Suyu (2003), and provides a 
good match with the observational data for M31. An A^- 
body realization of this model is done by GalactICS 
code written by Kuijken & Dubinski (1995). The satel- 
lite is assumed to be a Plummer sphere with initial mass 
10^ Mq (Model A), 5 X IO^Mq (Model B) and 10^° Mq 
(Model C). The initial scale radius is 1 kpc. Follow- 
ing Fardal et al. (2007), we adopt the initial position 
vector and velocity vector for the standard coordinates 
centered from M31 are (-34.75,19.37,-13.99) kpc and 
(67.34,-26.12,13.50) km s'^, respectively. The evolu- 
tion of the coUisionless system is followed up to 5 Gyr us- 
ing the parallel A^-body code Gadget-2 (Springel 2005) 
and AFD (Mori & Umemura 2006). We adopt a tree 
algorithm with a tolerance parameter of 6 = 0.5 and a 
softening length is 50 pc for all particles. The number 
of particles is 7, 351, 266 for the disk, 2, 581, 897 for the 



bulge, and 30, 750, 488 for the dark matter halo. For the 
satellite, we use 10^ particles for Model A, 5 x 10^ par- 
ticles for Model B, and 10^ particles for Model C. The 
total number of particles is « 40 million, and the mass 
of a particle is IO^Mq. 

Figure 1 shows the results for the time sequence of stel- 
lar dynamics for Model A. The upper panels in double 
columns illustrate the spatial distribution of stellar den- 
sity as a function of time from —1 to —0.25 Gyr (where 
Gyr is present-day). The lower panels describe the future 
evolution of the system from the present-day to 3 Gyr. 
In our time frame, —1 Gyr corresponds to the initial 
condition and the start of the simulation run. The first 
pericentric passage occurred about 0.8 Gyr ago. Figure 
lb shows that the satellite collides almost head-on with 
the bulge. Then, the distribution of satellite particles is 
distorted and is spread out significantly as seen in Fig. 
Ic. A large fraction of the satellite particles acquires a 
high velocity relative to the center of M31. This debris 
expands at great distance, remaining coUimated and it 
gives rise to the southern giant stream. The apocentric 
passage occurred ~ 0.7 Gyr ago and the second peri- 
centric passage is shown in Fig. Ic. As seen in Fig. 
Id, stellar particles that initially constituted the satellite 
start to form a clear shell structure after the second col- 
lision with the disk. Moreover, a double shell system is 
sharply defined in Fig. le. The system is composed of 
approximately constant curvature shells formed by phase 
wrapping (Hernquist & Quinn 1988). Then, these shells 
further expand and a multiple large scale-shell system is 
finally formed in the outer region and the dense core is 
formed in the inner region. The outermost large-scale 
shells in our simulation have a radius of > 50 kpc and 
these structures survive at least 4 Gyr from the present- 
day. 

Figure 2 shows the comparison between the map of 
RGB count density around M31 observed by Irwin et al. 
(2005) and the projected stellar density at the present- 
day in our simulation. The satellite is entirely disrupted, 
and the giant stream of debris arising from the tidal de- 
struction of the accreting satellite at the southern part 
of M31 is observed. The total mass of the stream given 
by the simulation is 1.4 x 10^ Mq. This is consistent 
with the estimated mass ^ 3 x IO^Mq from the obser- 
vations by Ibata et al (2001). The maximum length and 
width of the present-day model stream are ^ 150 kpc 
and ^ 50 kpc, respectively, in good agreement with the 
observations of Ibata et al. (2007). Furthermore, the 
simulation reproduces the butterfly-shaped shells in the 
northeast and the west part of M31. In the RGB count 
map of Irwin et al. (2005), the brightness of the north- 
east shell appears to exceed that of the west shell. The 
observed density contrast between the shells is 0.43 at 
(^,7?) = (2.0,0.7) and (-1.5,0.7) with the radius of 0.5 
deg for the standard coordinates centered from M31. On 
the other hand, the simulated density contrast is 0.25, 
which reasonably matches the observed faint features in 
M31. 

Figure 3 shows the disk thickness < z >^^^ as a func- 
tion of the radius following the interaction, for different 
satellite models. The solid line and the dashed line show 
the initial radial profile and the resultant radial profile 
after 5 Gyr without the infalling satellite, respectively. A 
concordance of two lines confirms that our fiducial model 
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Fig. 1. — Morphological evolution of the satellite, (a)-(f) The spatial distribution of stellar mass density as a function of the elapsed 
time, where 0.00 Gyr corresponds to the present. Each simulation box has a size of 160 kpc. The color corresponds to the logarithmic 
density 3.5 < logp (Mq kpc~ ) < 4.5. Elapsed time in Gyr is given at the lower-right corner in each panel: (a) — 1.00 Gyr, (b)— 0.75 Gyr, 
(c)-0.50 Gyr, (d)-0.25 Gyr, (e) 0.00 Gyr, (f) 1.00 Gyr, (g) 2.00 Gyr, and (h) 3.00 Gyr, respectively The formation of shells in the panel 
(e)-{h) infers that the distant arcs observed by Ibata et al. (2007) may correspond to the shells of an ancient radial merger. This figure is 
also available as an mpeg animation in the electronic edition of the Astrophysical Journal. 
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Fig. 2. — Gomparison between our simulation and the observed 
starcounts. Upper: The map of RGB count density around M31 
observed by Irwin et al. (2005). Lower: The projected stellar 
density of satellite particles for Model A. 



of M31 has perfect dynamical stabihty during the inte- 
gration time. There is no significant impact on the disk 
kinematics or disk thickness for Model A and B. There- 
fore, satellites less massive than IO^Mq have a negligible 
effect on the disk dynamics. But it is clear that the mas- 
sive satellite (Model C) more effectively heats the disk 
than do the less massive satellites. For the inner part of 
the disk, i? < 10 kpc, the scale height increases less than 
~ 10% of its initial value for Model A and B, while it 




R (kpc) 

Fig. 3. — Vertical extent of the disk at the present-day for 
different impact satellite masses. The solid line and the dashed 
line show the initial radial profile and the resultant radial profile 
after 5 Gyr without a collision, respectively. The dot-dashed line, 
the dotted line and the dot-dot-dot-dashed line correspond to the 
satellite mass of W^Mq (Model A), 5 X W^Mq (Model B) and 
10^'' Mq (Model C), respectively. We notice that our preferred 
Model A has a negligible effect on the disk. 



increases more than ^^ 50% for Model C. Thus, a mas- 
sive satellite is ruled out as the progenitor of the giant 
stream, and this result confirms the discussion in §2. 

4. DISCUSSION 

In this work, we follow the approach of earlier studies 
(Ibata et al. 2004; Font et al. 2006; Gechan et al. 2006; 
Fardal et al. 2006, 2007) in simulating the infall of a mas- 
sive satellite into M31. However, we add a live A^-body 
disk and bulge, and this permits us to ask under what 
conditions this collision might eject stars from these pop- 
ulations into the halo and quantifies the effect of these 
populations on the satellite. We also follow the evolution 
of this encounter for 4 Gyr beyond the present. As is the 
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Fig. 4. — Phase space distribution of the satellite particles. The 
color of each particle corresponds to the initial total energy with 
the unit of 10^^ erg, with more negative values having originally 
been deeper in the potential. Note that most of the particles that 
were originally deep in the potential of the satellite are found at 
present day in the East shell. 



case in Fardal et al. (2007), our simulations successfully 
reproduce the giant stream and (for a similar but much 
older collision) the apparent shell structures that are ob- 
served in the starcount maps. The large number of test 
particles in our satellite affords a good comparison with 
observations (Koch et al. 2007). We conclude that the 
event responsible for the giant stream most likely has not 
populated the halo with ejected disk stars; hence the in- 
termediate age metal rich populations observed in deep 
HST fields by Brown et al. (2006) are not disk stars 
ejected by the stream progenitor. An ancient radial col- 
lision involving a more massive progenitor (w IO^^Mq), 
or a collision with a different approach angle, might have 
been capable of ejecting disk stars into the halo. The 
present-day thickness of the disk following this collision 
constrains the progenitor mass to be < 5 x 10^ Mq. 

Our simulation gives additional predictions, such as 
the existence of high velocity stars in the central region 
of M31. Figure 4 shows that the tidal debris in the M31 
central bulge has velocity > 300 kms^^. The particles 
that were originally deepest in the potential (see color- 
coding in 10^^ erg units) are likely, at present, to lie in 
the East shell. Observations in this region might reveal 
additional clues about the nature of satellite, such as the 
central core and the metallicity gradients (cf. Mori et 



al. 1997, 1999). It is also noteworthy that > 3 Gyr 
after the collision, the debris more or less uniformly fills 
the halo. Simulations such as Font et al. (2007) may 
wish to consider the disk, bulge, halo, and satellite stellar 
populations and their effect on the long term evolution 
of streams; encounters with stars in the disk and bulge 
may disrupt the coherence of streams and scatter the 
stars in a manner analogous to disk and bulge shocking 
of globular clusters. 

A series of arc-like structures discovered by Ibata et 
al. (2007) are similar to the large-scale shell structures 
shown in Fig. If-lh. We conclude that these arcs may be 
the fossils of previous radial mergers several Gyr in the 
past (cf. Hcrnquist & Quinn 1988). Our model predicts 
that similar arcs should be found on the opposite side 
of the M31 disk, in the Northwest quadrant, if in fact 
they are the shells of ancient radial infall collisions. Our 
model also offers a natural explanation for stars with 
R > 100 kpc that are observed in the M31 halo near the 
systemic velocity of M31 (Gilbert et al. 2006; Koch et 
al. 2007). Some of these stars at i? « 165 kpc are found 
to have [Fe/H]=-1 by Koch et al. (2007). We propose 
that these stars are not a pressure supported halo but 
are rather the remnants of ancient radial collisions (note 
that in Fig. 1, the final extent of the collision is of order 
100 kpc). Future studies will examine collisions from 
different attack angles, stellar disks, and the fate of dust 
and gas. 
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